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Functional materials design normally focuses on structurally-ordered systems 
because disorder is considered detrimental to many important physical prop¬ 
erties. Here we challenge this paradigm by showing that particular types of 
strongly-correlated disorder can give rise to useful characteristics that are in¬ 
accessible to ordered states. A judicious combination of low-symmetry build¬ 
ing unit and high-symmetry topological template leads to aperiodic “ procrys¬ 
talline ” solids that harbour this type of topological disorder. We identify key 
classes of procrystalline states together with their characteristic diffraction be¬ 
haviour, and establish a variety of mappings onto known and target materials. 
Crucially, the strongly-correlated disorder we consider is associated with spe¬ 
cific sets of modulation periodicities distributed throughout the Brillouin zone. 
Lattice dynamical calculations reveal selective disorder-phonon coupling to 
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lattice vibrations characterised by these same periodicities. The principal ef¬ 
fect on the phonon spectrum is to bring about dispersion in energy rather than 
wave-vector, as in the poorly-understood “waterfall” effect observed in relaxor 
ferroelectrics. This property of procrystalline solids suggests a mechanism 
by which strongly-correlated topological disorder might allow new and use¬ 
ful functionalities, including independently-optimised thermal and electronic 
transport behaviour as required for high-performance thermoelectrics. 

The relationship between building-block geometry and bulk material structure is one of the 
cornerstones of structural science. By way of example, the solid phases of elemental Xe (1), 
Cgo (2) and human adenovirus (3) are all structurally related, not by virtue of any particular 
chemical similarity but because each of these phases reflects the same solution to the problem 
of packing weakly-interacting spherical objects in three-dimensional (3D) space. The reticular 
approach to understanding zeolite and metal-organic framework (MOF) topologies is related, 
because it links the geometries of molecule-like components (coordination polyhedra, molecu¬ 
lar linkers) to the 3D architectures formed by their assemblies (|?|[5j). The importance of struc¬ 
ture in determining the physical properties of solids is what then gives sense to the approach of 
developing new types of functional materials through informed design of their constituent build¬ 
ing blocks. It is precisely this type of “ground-up” approach that has recently been exploited 
in the rational design of e.g. solid oxide fuel-cell cathodes and room-temperature multiferroic 
candidates (6]|7j). 

To a large extent, the rich structural information accessible using crystallographic techniques 
has focussed effort on the design of crystalline materials. There are obvious functional advan¬ 
tages to the long-range periodicity characteristic of crystals, because it governs useful correlated 
properties—including the lattice dynamics, electronic states, and charge, orbital and magnetic 
order. Moreover crystal symmetry is central to mechanical properties such as piezoelectricity 
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and ferroelasticity, and is clearly pivotal in determining phase transition behaviour. While cer¬ 
tain building block geometries allow—or can even force (jS|{P|)—non-crystalline assemblies, the 
link to function is usually much less clear in these cases. Indeed the received wisdom is that 
disorder is something to be avoided, despite increasingly strong empirical evidence that links 
disordered states to advanced functionalities m- So the development of disorder-property 
relationships, and the eventual control over these properties through suitable building-block 
design have emerged as key challenges in the field. 

Here we develop an approach of intentionally designing functional disordered materials by 
focussing on systems in which structural disorder is extremely strongly correlated. The type 
of disorder we consider is similar to that found in ice, and our paper begins by developing a 
generalisation of ice-like states to arbitrary materials geometries. We proceed to establish a 
link between the geometry of structural building blocks and the propensity for specific types 
of strongly-correlated disorder in the resulting material assembly. This is the design element 
of our approach. Having made this connection, we develop a geometric basis for describing 
correlated disordered states. Our paper concludes by demonstrating how the correlated disorder 
deliberately engineered within one such material affects its lattice dynamics in a highly specific 
manner. This is the functional element of our approach because the effect we observe suggests, 
for example, a fundamentally new way of optimising thermoelectric response. 

Our starting point is the simple toy model of square ice, in which water molecules are 
arranged on a square lattice and oriented so as to satisfy sensible hydrogen-bonding rules 
[Fig. 1(a)]. There being no unique way of satisfying these rules the system is disordered, 
and—like the real-world examples of cubic ice (I c ) and its nano-confined variants (IT )—is char¬ 


acterised by a degenerate manifold of structural ground-states 02). An idea we will come to 
develop is that this propensity for disorder is encoded in the combination of the symmetry of the 
water molecule (i.e. the structural building block) and the lattice on which the water molecules 
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are arranged (here enforced by the directionality of the chemical interactions between building 
blocks). Any system that shares these geometric features will be characterised by the same con¬ 
figurational degeneracy. So, for example, replacing O-H.. .0 linkages by the M-C-N-M motif 
found in transition-metal cyanides gives a mapping that—in 3D—relates head-to-tail cyanide 
disorder in Cd(CN) 2 to water molecule orientations in cubic ice [Fig. 1(b)] ( 13). The ques¬ 
tion of 0/N ordering within square-grid-layers of transition-metal oxynitrides presents a related 
problem, which maps onto the square ice model following geometry inversion from one site to 


the next [Fig. 1(c)] (14). These examples involve compositional or orientational modulations 
of the square lattice, but the same ideas are well known to translate to a variety of modulation 
types, many of which are key to material function: e.g., displacive, electronic, charge density, 
spin density, orbital and spin orientation [Fig. 1(d)]. For ice-like disorder on the diamond lat¬ 
tice, these mappings are well established in the literature: hence the “Coulomb phases” (15) of 
charge (16), orbital (j77]), and spin (J7^) ices. 

In seeking to generalise ice-like states, we take our lead from the reticular chemistry ap¬ 
proach for generating network structures (79j 20). The key idea here is that lattice topologies 
can be considered in terms of the assembly of nodes and linkers [Fig. 2(a)]: the geometry of 
the node determines the possible topologies of the corresponding lattice. In this context, the 
various square-ice systems of Fig. 1 can be considered perturbations of the square lattice in 
which two adjacent linkers are distinguished amongst the four that meet at each square node. 
If we identify three linkers rather than two, then we generate a distinct family of disordered 
configurations that—from a reticular chemistry perspective—might be considered to arise from 
the linking of T-shaped building blocks [Fig. 2(b)]. As for the square ices, there are many 
possible realisations of this same state: one mapping is to Anderson’s resonance valence bond 


(RVB) description of singlet pair formation in cuprate superconductors (27); another is to so- 
called “domino” tilings of the plane d22j). There are in total just six cases to be considered for 
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perturbations of the square lattice. Two of these are trivial (distinguishing either four or zero 
linkers); two are related to one another (distinguishing one linker being the same as distinguish¬ 
ing three); and the two cases that remain distinguish different pairings of the four linkers, as 
shown in Fig. 2(c,d). We have already met the first of these cases in the guise of square ice 
[Fig. 1(c)]; the second case—in which the linkers distinguished are opposite one another—is 
ordered and results in symmetry breaking of the underlying square lattice. So there is a nontriv¬ 
ial relationship between perturbations of the node symmetry and the resulting configurational 
degeneracy. 

Ice-like configurations are not confined to perturbations of the square lattice. Equivalent 
states for the hexagonal, triangular, diamond, and cubic nets are enumerated in Fig. 2(e)- 
(s); those for the pyrochlore lattice are included as SI. The extent of disorder can be deduced 
from the corresponding diffraction patterns, which contain structured diffuse scattering in cases 
where there is strong correlated disorder [Fig. 2] (10 j). What emerges is that a substantive frac¬ 
tion of these systems admit large configurational entropies, with a complex relationship between 
node geometry and extent of correlated disorder. So as to provide insight into this relationship, 
we generalise Pauling’s approximation for the configurational entropy of ice (12 [74]): 
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Here d represents the underlying node connectivity and n corresponds to the number of symmetry- 
equivalent node perturbations of a given type. The value of n can often be determined by in¬ 
spection, but it is given more rigorously by the ratio of the orders of the point groups of the 
parent and perturbed node geometries; for example, n = \D 4 h \/\C 2 v\ = 4 for square ice. We 
call p = n/2' l l 2 the Pauling number, with the significance that maximising p maximises the 
propensity for disorder. In this way one expects low-symmetry perturbations of high-symmetry 
lattices to lead to states of the greatest configurational entropy, a qualitative relationship that is 
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borne out in practice [Fig. 2] .From a materials design perspective, what we are saying is that 
building block geometry and the arrangement of the interactions between building blocks can 
together encode for specific types of correlated structural disorder. 

Common to many of the various configurations of Fig. 2 is the absence of translational 
periodicity characteristic of the crystalline state. For a given configuration, every node experi¬ 
ences the same local environment—and hence it is not meaningful to think of these structures 
as defective in the vernacular sense. Yet there is no unit cell and space group symmetry that 
properly describes the topological connectivity. We proceed to argue that these systems should 
not be considered as crystals, but form a separate class of aperiodic solid with its own char¬ 
acteristics. We will use the term procrystalline to describe this state and to emphasise that 
conventional crystals might be seen as a special case of the definitions that follow. The procrys¬ 
talline state is a dense, overlapping packing of identical fundamental structural units (we use 
the term “ neighbourhoods ”), which are positioned periodically but orientationally permuted as 
permitted by the point symmetry of the neighbourhood geometry. For magnetic systems, these 
permutations may involve time reversal operations as realised in e.g. the Ising spin ices (23). 
In simple cases the neighbourhood corresponds to the Dirichlet-Voronoi cell of the underlying 
lattice augmented to include neighbouring, correlated lattice points [Fig. 2(b)] (see SI for fur¬ 
ther discussion). Whereas crystals correspond to the special case in which the neighbourhood 
orientations are themselves periodic, the more general procrystalline state allows for discrete 
orientational disorder (c/. the continuous orientational degrees of freedom in e.g. plastic crys¬ 
tals and superionics). Any such disorder will always be correlated since neighbourhoods over¬ 
lap. Because the underlying neighbourhood lattice is periodic, all procrystals admit a Bragg 
diffraction pattern and have a well-defined reciprocal lattice. This diffraction pattern can be 
analysed using conventional crystallographic approaches but doing so yields a structural model 
in which neighbourhoods are averaged over their different possible orientations and all infor- 
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mation regarding orientational correlation is lost; for example, the states represented by panels 
(r) and (s) in Fig. 2 share identical Bragg diffraction patterns in spite of their distinct local sym¬ 
metries. Like crystals, procrystalline phases are characterised by macroscopic point symmetry 
that can be as high as that of the neighbourhood lattice. Yet, unlike crystals, they can support 
a complete absence of any point or translational symmetry at the microscopic level. It is the 
existence of a periodic 3D reciprocal lattice that distinguishes procrystals from incommensurate 
and quasicrystalline phases, and which also guarantees a well-defined Brillouin zone (BZ) and 
Bloch-wave description of phonon and electronic states. 

The structures of a number of well- and lesser-known materials can be thought of in pre¬ 
cisely these terms. This is true by construction for any phase with ice-like disorder; in addition 
to the various systems described above, the family of ferroelectric phases related to KH 2 P0 4 
is an obvious additional example (j 24). Similarly well established are the statistical mechani¬ 
cal models of RVB (J25|) and loop d23j) states, which correspond to procrystalline lattices with, 


respectively, one and two linkers distinguished for each node. Hence physical realisations of 
either class also fall under our definition (e.g. TaS 2 f27j) and SrTa0 2 N (28)). A less obvi¬ 
ous example is the assembly of p-tcrphcnyl^.SA'.S'-tctracarboxylic acid (TPTC) molecules on 
pyrolytic graphite to form a hydrogen-bonded network related to the procrystalline lattice il¬ 


lustrated in Fig. 2(e) [Fig. 3(a,b)] (29). This arrangement maps onto the so-called “rhombus” 
or “lozenge” tiling, which in turn corresponds at once to both the Ising triangular antiferro- 
magnet and the RVB description of 7r-bonding in graphene (j30 31 ~32\ . These equivalences 
are straightforwardly seen in reciprocal space: Fourier transform of the scanning tunnelling 
microscopy image of Fig. 3(b) reveals the same distribution of diffuse scattering and “pinch 
point” features expected from our simple geometric model [Fig. 3(c,d)]. Indeed we expect the 
link between characteristic diffuse scattering patterns and particular procrystalline states to aid 


in diagnosing and understanding a range of problems of correlated disorder (10 33). For pow 


7 








der diffraction measurements, often the only signature of this diffuse scattering is the presence 
of /^'/-dependent anisotropic peak shape broadening. This is the case, for example, in the scat¬ 
tering patterns of Pd(CN) 2 and Pt(CN) 2 ; a procrystalline structural model based on connected 
square-planar [M(C/N) 4 ] units provides the first convincing description of their diffraction be¬ 


haviour [Fig. 3(e,f)] (34). In other cases, procrystalline states (if not necessarily recognised as 
such) have been inferred from the combination of (i) a disordered average structure, with (ii) 
clear signatures of local distortions that can persist only if suitably correlated. Examples include 


Jahn Teller distortions in the high-temperature orbital-disordered phase of LaMn0 3 (35 36) and 


the high-pressure amorphous phase of ZrW 2 0 8 p7j). So there is good evidence that a variety of 
procrystalline phases do exist, even if their structures are difficult to interpret using established 
crystallographic approaches. 

But what of the link between correlated structural disorder and function? In principle, the 
existence of a well-defined BZ allows coupling between the structural modulations that char¬ 
acterise the procrystalline state and other physical properties that depend on periodicity— e.g. 


the lattice dynamics and electronic band structure (38). We tested for coupling of this type 
using as our example a two-dimensional oxynitride lattice [Fig. 1(c)]. The idea was to set up 
a simple harmonic lattice dynamical model in which we assigned different equilibrium val¬ 
ues and stiffnesses to M-0 and M-N bond lengths, and again to O-M-O, N-M-O, and N- 
M-N bond angles, and then to determine the extent to which correlated compositional order 
affected the phonon spectrum. Conventional lattice dynamical calculations are designed for pe¬ 
riodic (crystalline) structures, and so we made use of a supercell Monte Carlo lattice dynamical 


(MCLD) approach (39 40) in order to treat compositional disorder explicitly (see SI for further 
discussion). We first benchmarked our calculations by determining the “mean-field” phonon 
dispersion expected for an average of the different force constant values; we found essentially 
perfect agreement between our MCLD calculations and those obtained using the GULP pro- 
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gram (41) [Fig. 4(a)]. For random distributions of equal numbers of O and N atoms, the basic 
phonon structure was similar to the mean-field case, with a slight redistribution of energies 
among phonon branches [Fig. 4(b)]. This result is in agreement with ab initio molecular dy¬ 
namics studies of the related problem of the lattice dynamics of configurational glasses (42 j). 
By contrast, for the procrystalline arrangement there was a dramatic dispersion in energy of 
the optic branches that we observed only for the specific subset of wave-vectors of the form 
(4, £)*—precisely the family of modulation periodicities characterising the procrystalline state 
in this case [Fig. 4(c,d)]. The qualitative similarity to the “waterfall” phonons observed in 
thermoelectrics and relaxor ferroelectrics is striking, and suggests a plausible origin for the 


phenomenon in those systems (43 44). 

So this is our key result: strongly-correlated structural disorder allows selective control over 
physical properties that depend on periodicity. The implication for systems where thermal and 
electronic conductivities are mediated by, respectively, phonons and electronic states localised 
in different regions of the BZ is that disorder-phonon coupling offers a means of selectively 
reducing thermal conductivity (inversely proportional to phonon bandwidth) without affecting 
charge transport behaviour. This is an attractive design strategy for developing next-generation 
thermoelectrics, and one that contrasts with the use of “rattlers” which are indiscriminant in 
their k-space coupling ( |45] ). 

Our reticular chemistry methodology suggests a number of synthetic routes for realising 
new classes of functional procrystalline solids. Metal-organic frameworks (MOFs) are an obvi¬ 
ous platform, given (i) they offer the requisite control over building unit geometry and (ii) their 
energetics tend to be dominated by local interactions (46 ?7j). While there is reduced scope for 
coupling between structural disorder and electronic behaviour in these systems, porosity perco¬ 
lation will certainly be affected by disorder and may in turn govern sorption, mechanical and 


ion storage properties (48 49). In more conventional inorganics, local symmetry lowering can 
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be achieved by covalency effects (as in mixed-anion perovskites) or by first- or second-order 
Jahn Teller distortion (as in the chalcogenide thermoelectrics). Moreover, because our analysis 
is essentially geometric in nature, there is clear scope to extend these concepts to magnetic or 
electronic states, or indeed to the macroscopic scale. The recent demonstration that disordered 
metamaterials can show strong structural coupling to light scattering processes is an example 
relevant to the generation of modern photonics (50). Thinking beyond the ground-state proper¬ 
ties of procrystals, we anticipate the existence of novel collective and hidden degrees of freedom 


that promise a rich physics of their own; e.g. , topological excitations (51) and/or “hidden order” 


transitions between distinct local symmetries (52). 
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Fig. 1. 

Correlated disorder in square ice analogues, (a) A configurational fragment of square ice. 
The water molecules are arranged on a square lattice and are oriented so as to satisfy local 
hydrogen-bonding rules: each molecule accepts two hydrogen bonds and donates two hydro¬ 
gen bonds. There is no unique solution to satisfying these local constraints, and so the square 
ice state is configurationally disordered, (b) A square-planar transition-metal cyanide config¬ 
uration that maps onto the ice-like state in (a). Here each metal cation (coloured green) is 
coordinated by two nitrogen atoms (blue) and two carbon atoms (grey) such that each N atom is 


opposite to an C atom, (c) Transition metal oxynitrides adopt a related structure U4), in which 
square-grid layers consist of metal cations coordinated by two nitrogen atoms (blue) and two 
oxygen atoms (red). The topological equivalence to the square ice configuration can be seen 
by alternately considering O-M-O and N-M-N orientations for neighbouring metal centres 
(shaded regions) (2<5J). (d) Displacive modulation of a square M0 2 lattice by in-plane [M0 4 ] ro¬ 
tations gives configurations that again map onto the square ice state. The correspondence relates 
displacements above and below the plane with, respectively, O and N atoms of the oxynitride 
configuration shown in (c). 
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Fig. 2. 

Reticular design approach to generating procrystalline networks, (a) High-symmetry build¬ 
ing blocks connect to form familiar two- and three-dimensional networks: square (indigo), 
hexagonal (blue), triangular (green), diamondoid (orange), and cubic (red), (b)-(s) Distin¬ 
guishing different possible subsets of linkers for these high-symmetry lattices gives a variety of 
ordered and disordered states, which are grouped here according to the parent lattice. For each 
panel, the perturbed node geometry is shown in the top-left corner, followed immediately below 
by a representation of the Pauling number p, a qualitative indicator of the propensity for disor¬ 
der. A representative network configuration is shown as the main image, with nodes coloured 
according to their orientation. One suitable projection of the corresponding X-ray diffraction 
pattern is given in the bottom-left corner (see SI for further details). For the configuration shown 
in panel (b), two overlapping neighbourhoods are outlined in black. The Dirichlet-Voronoi cell 
of the neighbourhood lattice is shown in red; the neighbourhood itself is generated by augment¬ 
ing this cell to include connected latticed points. The asterisk in panel (1) indicates that a single 
enantiomer of the node geometry is used (c/. panel (k)); this node is chiral when constrained to 
lie in two dimensions. Further discussion, including extension to the pyrochlore lattice, is given 
as SI. 
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Fig. 3. 

Representative physical realisations of procrystalline states, (a) Molecules of p-tcrphcnyl- 
AAA'A'-tctracarboxylic acid (TPTC) self assemble on pyrolytic graphite to form a hexagonal 


procrystalline network related to that represented in Fig. 2(e) (29). One possible arrangement 
is shown here (left), together with the corresponding topological representation of the type 
used in Fig. 2 (centre); molecule orientations map onto the “missing” linkers of the topological 
representation. Shown in background is the corresponding rhombus tiling, which shares the 
same configurational degeneracy. There is a further one-to-one correspondence to ground states 
of the triangular Ising antiferromagnet (right): opposite vertices of each rhombus are decorated 
with pairs of “spin-up” (white circles) and “spin-down” (black circles) states. The constraint 
that each triplet of neighbouring spins (red triangle) contains at least one spin-up and one spin- 
down state perturbs the node symmetry of the underlying hexagonal lattice in the same way as 
does the pattern of hydrogen- and covalent-bonding in the molecular assembly, (b) A scanning 


tunnelling microscopy image of the corresponding experimental state [Adapted from Ref. (29)]. 

(c) The Fourier transform of the image in (b), shown regions of structured diffuse scattering. 

(d) The strikingly similar scattering pattern calculated for the hexagonal procrystalline network 
of Fig. 2(d). (d) A procrystalline model for the structure of Pd(CN) 2 and Pt(CN) 2 based on the 
state represented in Fig. 2(p): the two possible orientations of square-planar M(C/N) 4 nodes are 
shown in indigo and gold, (f) Comparison of the Rietveld fit for this structural model of Pd(CN) 2 
(red lines) and the experimental X-ray powder diffraction data of Ref. p9) (black points, A = 
1.54 A); tick marks indicate the positions of parent Bragg reflections and the difference (data- 
fit) is shown in blue (see SI for further discussion). 
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Fig. 4. 

Selective disorder-phonon coupling in a procrystalline network, (a) Mean-field phonon 
dispersion curves for a two-dimensional oxynitride lattice, determined using MCLD (coloured 
bands) and GULP (black lines) approaches. The width of the coloured bands reflects the stan¬ 
dard uncertainty of the MCLD phonon frequencies, which is non-zero as a result of the finite 
configurational ensemble size used for phonon calculations (see SI for further discussion), (b) 
The phonon dispersion determined using MCLD for configurations in which equal numbers of 
O and N atoms are distributed randomly across the O/N sites of the same oxynitride lattice. 
There is a slight redistribution of phonon frequencies relative to (a) but the standard uncertain¬ 
ties remain essentially unchanged, (c) A third set of phonon dispersion curves, again determined 
using MCLD but for configurations in which O and N atoms have been distributed according to 
the correlated disorder (procrystalline) model of Fig. 2(c). There is now substantive “waterfall”- 
like phonon broadening along the X-M direction which is an order of magnitude larger than 
the intrinsic finite-ensemble effects, (d) The calculated single crystal scattering pattern for the 
procrystalline model of Fig. 2(c), demonstrating that this particular state is characterised by 
modulations that—like the phonon broadening—are localised along the X-M branch of the BZ. 
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